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Abstract 

! Although dark matter makes up 80% of the gravitational mass of our Galaxy, its composition 

is not known. One hypothesis is that dark matter consists of massive particles called WIMPs. 
WIMPs are expected to accumulate and coannihilate in the cores of stars, but the only signature 
of this accumulation has been thought to be hard- to-observe high-energy neutrinos. Here we 

& : 

propose an entirely new observable signature. WIMP coannihilations in the core of a very low- 



mass star, brown dwarf, or planetary-mass object should alter the star's chemical composition via 



■ spallation reactions. Very close to the Galactic center, these stars may acquire extremely high 

lithium, beryllium, and boron abundances, even for models with otherwise- undetectable wimp- 
co 

nucleon cross sections. These abundances should be measurable in certain stellar systems and 



phenomena. 
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I. INTRODUCTION 



In this era of high precision cosmology, the mass-energy density of the Universe is well- 
constrained. From a combination of measurements of the microwave background, large- 
scale strucure, rotation curves, gravitational lensing, and galaxy cluster temperatures, the 
Universe is known to contain about 4% ordinary matter and 24% dark matter. Dark matter 
is known to move and gravitate just like massive baryonic matter, but without any particle- 
particle interactions at the electromagnetic or strong-force scales. The Standard Model of 
particle physics contains no massive non-interacting particle like this, but some extended 
models do; such particles are called WIMPs. Supersymmetry (SUSY) naturally includes 
a WIMP, the neutralino (xo), which is an excellent dark-matter candidate; Kaluza-Klein 
bosons, which arise in theories with large or warped extra dimensions, are another example. 
In this work we refer to WIMPs as Xo, but our results are not model-specific (except as 
discussed below). 

In order for WIMPs to be produced in the Big Bang, there must be a sizeable cross section 
for reactions like xx — > XoXo ; where x is some Standard Model particle. This implies that 
the reverse reaction, XoXo —> will occur today at a rate proportional to the square of the 
WIMP density. The WIMP-nucleus elastic scattering cross section, XqN z - > Xo-^z is n °t 
constrained by astrophysics, and may be extremely small. Terrestrial experiments allow us 
to constrain this cross section to be <7o/A 2 <10~ 43 cm 2 for spin-independent WIMP-nucleus 
scattering, and o" p <10~ 36 cm 2 for spin-dependent WIMP-proton scattering. 

Any nonzero cross section for Xo^z ~~ > XqN z allows a Galactic Xo particle to scatter off 
nuclei in stars, to enter bound orbits which cross the star repeatedly, and thus to thermalize 
and be captured inside the star. As the xo number density in the star's core increases, the 
coannihilation rate will increase until it equilibrates with the capture rate. Recent work l| 
shows that, in regions near the Galactic center's SuperMassive Black Hole (SMBH) where 
the dark-matter density is highest, the capture- and- annihilation process may provide enough 
thermal power to change a star's luminosity and evolution. We show that these coannihi- 
lations can change stellar chemical composition in unambiguous ways. Most importantly, 
the high-energy annihilation products can convert carbon and oxygen into 3 He, Li, Be, and 
B via spallation. We restrict our analysis to stars with M < 0.1 M Q (encompassing the 
lowest-mass main-sequence stars, brown dwarfs, and planetary-mass objects) which we will 
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FIG. 1: Spallation nuclei produced per primary GeV, based on simulations of b-, W-, or t-jets 
incident on degenerate SSO core matter. Yields for mixed solar material are consistently a factor 
of ~ 100 lower. 

refer to loosely as substellar objects (SSOs). 

II. SPALLATION REACTIONS AND WIMP CAPTURE 



In order to quantify the spallation yields, we perform a numerical simulation of WIMP 
coannihilation products colliding with stellar core matter. Many SUSY WIMP models sug- 
gest that XoXo - * xx usually proceeds, via prompt heavy quark or W boson pairs, into 
hadron jets containing 50-100 particles. Using the Pythia6 event generator |2j, we create 
WIMPs at rest and annihilate them to the heaviest allowed of bb, W + W~, or ti\l^ . We 
use the Geant4 Monte Carlo package ^] to propagate the final states through two relevant 
compositions: (i) a uniform solar composition [Z =1-26); and, (ii) the approximate core 
composition of a gravitationally-differentiated star (solar composition with Z =6-26) 20]. 
The spallation yield per unit xo mass is shown in Figure [1] for the core material. Note that 
the yields for a fully mixed composition are a factor of 100 lower. 

Next, we derive the WIMP capture rate for a range of example stars near the Galactic cen- 



We assume that the Galactic 



ter, using the analytical formula (equation 2.27) derived in 
halo has an adiabatic central "spike" with a static density profile of p(r) = lO 5 M /pcxr~ 2 . 
Many dark-matter halo profiles are assumed to have a spike |5J] as long as it has not been 
disrupted by a galaxy merger; our assumption is consistent with, but not required by, the 
observations. We have not accounted for the evolution of the spike, whose density would 
have been higher earlier in its formation history [6]. (To adapt our results to other halo/ spike 
models, the results may be scaled linearly with WIMP density.) To simplify the calculations, 
we treat the stars as n = 1.5 polytropes in circular Galactic orbits. There are two distinct 
possibilities for the WIMP elastic scattering cross section: (i) spin-dependent; and, (ii) spin- 
independent. Spin-independent capture on SSOs is very inefficient, since they consist mostly 
of light elements while the cross section scales as the nuclear mass squared. However, the 
spin-dependent interaction may permit efficient scattering off of spin- 1/2 protons (hydro- 
gen), allowing SSOs to be significant dark-matter burners. The capture rate depends on the 
SSO's mass (M*), WIMP mass (M Xo ), the distance from the Galactic center (R) and the 
WIMP-proton cross section. As long as o v < 10~ 36 cm 2 , this rate simplifies to 

r _ 1Q33 ' 6 r °P ) l( r ) -l/2 ( M * ] 3 ( R * ] ~2 ( M XQ n-2 m 

s l 10- 36 cm 2j l pc M Q } { R Q } l 100GeV J 1 1 

The capture rate is independent of a„ for o v > 10~ 36 cm 2 . 



III. SPALLATION OF SUBSTELLAR OBJECT CORES 

If the WIMPs are in thermal equilibrium with the star, they will be collect near the star's 
core. The creation of spallation products and the production of energy will also be largely 
confined to the core region. This has important implications both for the thermal[l] and 
the chemical structure of these stars. 

In a uniformly-mixed SSO, spallation would occur mainly on H and He, with B produc- 
tion on C and O being suppressed. However, many SSOs have likely experienced chemical 
differentiation wherein most of their metals would have gravitationally settled into a central 
coreQ; this core is large enough to enclose the entire WIMP annihilation region. In these 
differentiated stars, the spallation target will be mostly C and O (producing abundant Li, 
Be, and B), with traces of elements up to Fe (producing rare elements like Sc, V, and ra- 
dioactive 26 Al and 40 K). Nuclear re-burning of the spallation products is negligible for stars 
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TABLE I: Seven examples of Galactic-center SSOs. We compute the spallation-product contents 
for an arbitrary dark matter model with m Xo = 100 GeV and spin-dependent a p = 10 -38 . We 
calculate the thermal power and capture/burning rate at equilibrium. We also show the average 
boron abundance of a fully mixed SSO after 10 Gyr, first under the assumption that the SSO is 
gravitationally settled and second under the assumption that it is fully mixed. 

with M < O.1M . Stars with 0.07 < M/M < 0.1 can burn Li and Be but not B, and 
stars with M > 0.1M Q can burn B|8(. We note that, unlike a thermal signature, a chemical 
signature may persist indefinitely for a star which has been ejected from the Galactic center 
and has ceased burning WIMPs. 

For concreteness, we define a set of representative substellar objects and calculate their 
expected WIMP annihilation rates and spallation accumulation rates for a WIMP mass of 
M Xo = 100 GeV and a spin-dependent WIMP-proton cross section of <To = 10 -38 cm 2 . The 
results are shown in Table [B 

IV. OBSERVATIONAL PROSPECTS AND CONSTRAINTS 

Low-background nuclear physics experiments like CDMS [s| typically plot their dark- 
matter discovery potential as a function of a v and M xo . We do the same for spallation 
production in SSOs in Figure [2j The colored contours show where the example stars from 
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Table d can accumulate a net boron atomic abundance of > 10~ 8 in 10 Gyr. Our results 
demonstrate that extremely low-mass stars, brown dwarfs, and planetary mass objects can 
produce anomalous abundances of Li, Be, and B if: (1) these objects reside within the inner 
parsec of the Galactic center; (2) the dark matter consists of WIMPs with properties as 
described; and, (3) the WIMP density distribution has a spike at the Galactic center. Very 
few other assumptions are required, except that WIMP coannihilation proceeds to hadrons. 
The discovery or exclusion of anomalous B levels in stars like these could provide a discovery 
channel for, or powerful constraints on, various dark matter models that might otherwise 
be inaccessible. Moreover, any 'positive identifications' that may determined from other 
experiments will need to be checked by as many independent tests as possible. 

Unfortunately, the measurement of a detailed spectrum of even the brightest of these 
objects will be extremely challenging for any existing or foreseen telescope, due to the large 
distance (8 kpc), interstellar extinction, and crowding. However, some special cicumstances 
may allow us to observe the spallation byproducts of such objects in the future. For exam- 
ple, three-body encounters may be reasonably effiicient at ejecting low-mass objects from the 
deep gravitational potential well near the Galactic center and thus, on a gigayear timescale, 
seed the local environment with observable candidates. The detection of anomalous abun- 
dances of light elements such as boron in the atmospheres of isolated, very low-mass objects 
whose interiors are fully convective would be an important indicator of WIMP-induced spal- 
lation. 

Interacting binaries containing compact accretors, such as Cataclysmic Variables and 
Low-Mass X-Ray Binaries, can also expose the chemical constituents of the interiors of their 
low-mass companions (i.e., the donor stars). Thus our candidate objects may serve as the 
donors in these systems (e.g., CV SDSS 103533.03+055158.4 [h]]), and the accretion disks 
may make their composition visible in X-rays or the UV. Moreover, some of these systems 
can cannibalize their companions to such a degree that only the cores of the donors persist 
(e.g., the black widow pulsar or other ultracompact binaries). These types of binaries allow 
us to potentially observe the spallation products of very low-mass objects whose interiors 
had been chemically differentiated and whose core material could not have been otherwise 
dredged up to the surface. Moreover, if spallation products from the core have not mixed 
with the whole star, the core boron concentration may be orders of magnitude higher than the 
average concentrations discussed above. Finally we note that the total disruption of a very 
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low-mass object would give us the possibility of observing the core's chemical composition. 
In fact, one of our candidate objects may be disrupted by, and accrete onto, the SMBH 
itself (perhaps as in RX J 1242- 11 U|). The resulting burst may show evidence of spallation 
products. Observation of the boron abundances of Galactic Center SSOs or SSO-related 
phenomenon, although difficult, should be a high priority for future UV missions. 
The authors wish to thank Saul Rappaport, Adam Burgasser, and Gray Rybka. 
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FIG. 2: Discovery potential for dark matter in observations of spallation products. The graph 
shows the spin-dependent xo-P r °t° n elastic scattering cross section versus the WIMP mass. The 
colored solid contours show where the example stars defined in Table U would achieve a total B/H 
ratio > 10~ 8 (quantities in the legend are the mass in M and Galactic radius in pc). Examples 1 
and 2 assume complete mixing, all others assume full gravitational settling. Dashed (dotted) red 
lines are existing (proposed) experimental limits; top to bottom, they are: CDMsJjJ, KIMS [3], 



COUPpQ, NAIAD 1 I] (100 kg-yr). and Supcr-Kannokandc 
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